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Abstract 



Models for species-area relationships up to now have focused on the mean richness as a 
function of area. We present MFplp2, a self-similar multifractal. It explicitly models both 
trend and variation in richness as a fun ction of area , and is a generalisation of the model of 
scaling of mean species richness due to 



Harte et al 



( 1999 ). The construction is based on a 



cascade of bisections of a rectangle. The two parameters of the model are pi, the proportion of 
species that occur in the richer half, and p2, the proportion of species that occur in the poorer 
half. Equivalent parameterisations are a = (pi +p2)/2 and b = pl/p2. These parameters 
are interpreted as follows: a gives the scaling of mean density, b gives the scaling of spatial 
variability. Several properties of MFplp2 are derived, a generalisation is noted and some 
applications are suggested. 



INTRODUCTION 



Ecologists study organismal interaction on a wide range of scales. From the centimetres of a puddle 
to the megametres of the planet, the characteristic length of the system spans some nine orders 
of magnitude. The mere description of pattern over such a range of scales is a major challenge. It 
becomes formidable if in addition one requires that the description be simple. 

But species-area relationships (SARs) achieve exactly that. Let S denote the number of species 
found in a survey of a spatial domain with area A. Then one may hope that S increases with 
A in a regular way. I ndeed, und e r the interpretation that S{A) denotes the average richness of 



domains with area A 
set of islands, while 



Gleason 



Arrhenius ( 19211 ) found that the power law S — cA^ held for plants on a 



1922) found that the semi-log law S = a + b \og{A) held for plants 



in patches of a single vegetation type. Both of these laws have the form 

S^f{A), 

where the function / is explicitly given in a short formula with few parameters. 



(1) 



Since t hen many othe r likewise simple forms have been found to fit data on S versus A. The re- 
view bv lTjorvd (|2003l ) lists 14 such functions in its two main categories and mention s several more 



The m ost commonly used remain the power law and the semi-log law: for example, 



Drakare et 



(|2006t ) found 796 and 506 cases, respectively. Such a wealth of results demand interpretation, for 



which one requires a theoretical framework. Usually, this is along the lines elegantly presented 
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by iMavl ([l975|). The starting point is to assume known the species abundance distribution (SAD) 
and the total area (i.e. the largest area over which the SAR is to be applied). Next, it is assumed 
that the density of individuals is uniform throughout the area and that they are independently 
distributed. The expected number of species in a given area A then follows from a probabilis- 
tic calculation. In other words, this approach explains SARs from SADs. Normally the latter 
are derived from ecological processes using a second theoretical framework. One is thereby en- 
abled to interpret differences in the parameters or even in the form of equation [T] in terms of 
underl ying e cological processes. A high point of this approach is Rosenzweig's magisterial syn- 
thesis (|l995f) . which uses power law SARs interpreted in terms of island biogeography. Other 



theories for which this approach w orks include metapopulation d ynamics ( Hanski fc Gvllenbers 



1997 



Ovaskainen fc Hanski 



Chave . 



2004 ) 



2003| l and Hubbell's neutral theory (jVolkov et al 



2003 



Belli, 



200C 



A completel y diffe rent modelling paradigm is due to Harte and colleagues (jHarte fc Kinzig 



Harte et al 



1997 



19991 ). They invert May's sequence, and derive an SAD from an SAR. In their 
models, the spatial domain is a rectangle with area Aq, and it is repeatedly refined by a sequence 
of bisections into increasingly smaller rectangles. The SAR is applied at each bisection, thus 
specifying the number of species retained in each half. In their papers, Harte and colleagues use 
the power law S — cA^ , and this implies that the proportion of species retained in each half is 



a = 2- 



By further assuming that the bisection cascade ends with 2™ recta ngles of size Ao/2 



and that each such rectang le contains one individual, an SAD can be derived 



Martin fc Goldenfeldl . 



2002). 



Harte et al 



1999 



It should be noted that the SAD derived from S ~ cA^ in the Hartean scheme is not the SAD from 
which May derives S = cA^. The reason appears to be that in May's derivation each individual 
is indpendently placed, so that spatial autocorrelation is zero by assumption, while in Hartean 
models individual positions within species are correlated, in that individuals of a rare species are 
very likely to occur near each other, and hence spatial autocorrelation within species cannot be 
zero. 



The Hartean approach is attractive for several reasons. Firstly it directly gives species richness 
as a density which varies according to scale. Put differently, it provides a simple way to simulate 
species richness across many scales. Secondly, it is invariant with respect to scale. That is, 
the relationship between two patches depends only on their relative areas, not on their absolute 
areas. Thirdly, it is extremely simple. For example, it uses only one parameter, whereas May 
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uses two (actually three, but any one parameter can be derived from the other two). As a purely 
descriptive tool it has much to recommend it. However, its very simplicity counts against it. 
It has been pointed out that the bisection cascade creates a very limited set of area s, leading 



to th e claim that the model does not apply to areas not included in this small set (jMadduxl . 



2004L on which more below). Also, and perhaps more importantly, it is explicitly a model of 
average richness. This limitation is of course shared by all models in the form of equation [TJ 
Nevertheless, variabili ty in richness is a n obvi ous fact and indeed one of the main o bjects of 



study in biodiversity 



2005 



Pautasso , 



Connor &: McCov 



20011 ) and is known to be scale-dependen t 



20071 ). It also varies in time ( Gonzalez 



200C; 



Adler fc Lauenroth 



Rahbek 



200, 



description of species richness with appropriate scaling properties for the average richness as well 
as the variation in richness is likely to be very useful. 



In this paper we present such a model: MFplp2. It is a simple generalisation of the Hartean 
approach, where instead of assuming that both halves of a bisection retain the proportion a, 
we assume instead that they differ, with the richer half retaining the proportion pi and the 
poorer half retaining p2 of the species in the bisected rectangle. As a result, the rectangles with 
given area shows a range of richnesses. In other words, MFplp2 gives the SAR as a set-valued 
function: 5 at a given A takes a set of values. We give an exact expression for the frequency 
distribution of these values, and derive a number of consequences: self-similarity, a simple graphical 
characterisation, numerous power-law scalings, and the prediction that richness and endemism are 
positively correlated. 



Madduxl . 



2004 )? Do they not also 



What about objections to the original Hartean model (e.g. 
hold in this case? Yes, but we do not believe they invalidate the use of either MFplp2 or its 
Hartean parent. Firstly, these models are descriptive, not causal. Their value lies in their ability to 
summarise complex reality in a few parameters. They need merely do that with sufficient accuracy. 
In this respect, they resemble null models and neutral models. In spatial ecology at present such 
descriptions are indispensable, even when not very accurate. Without them, we have only the fully 
detailed maps of distributions and richness. Secondly, the limitations of the constructions used to 
explain these models can in many ways be mitigated by interpretative strategies and extensions. 
We explore both of these points more fully in the discussion. 
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Figure 1: Constructing MFplp2: stages fc = 0, 1, 2 and 3 of the multiplicative cascade. At each 
stage, there are 2^^ bisections, pi is the fraction of richness inherited by the richer half, p2 is the 
fraction of richness inherited by the poorer half. Position of richer half is randomly chosen among 
the two possibilities at each bisection. 



MFplp2: A TWO-PARAMETER MULTIFRACTAL MODEL 



This exposition refers to Figure [TJ Consider a rectangle with area Ao that contains 5*0 species. 
Bisect the rectangle into two congruent rectangles of area Ai — Ao/2, then bisect each of these 
with a line at right angles to the previous bisector to reach four rectangles of area A2 — Aq/A. 
Continue this process, reaching 2'^ rectangles of area Ak — Ao/2'' after k bisections. 

As each parent rectangle is split in two, spatial heterogeneity in species richness implies that 
its offspring inherit unequal amounts of species. Denote by pi the fraction of species inherited 
by the richer half and by p2 the fraction inherited by the poorer half. For example, after the 
bisection of ^0 the richer half contains piSo species, and the poorer half contains p2'S'o species. 
Since every species in the larger rectangle must occur in at least one of the smaller rectangles, we 
have pi + P2 > 1- At each stage, the position of the richer half is randomly assigned. Typical 
realisations are shown in Figure [2l 
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Figure 2: Realisations of MFplp2 with pi = 1, p2 = 0.6, for fc = 2, fc = 3, fc = 6 and k = 10. 
Density scale is the same in all four cases: the maximum richness is pf = 1, which appears as black 
areas (squares/rectangles) in each of the four examples, and the minimum richness is pl° w 0.006, 
which appears only in the A; = 10 example. 
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Notice that pi and p2 remain constant throughout the construction. After k bisections the het- 
erogeneity is captured in the binomial {pi + P2)''' — X]j=o "'?'2- Species richness takes the 
{k + 1) values p\^-' pi^So for j = to fc, and each level of richness occurs in (^) of the 2'^ rectangles 
with area If one denotes by Pk{S — n) the proportion of rectangles with area 2~^Aq that 
contain n species, then we can restate this result as the species richness distribution 



Pk{S = Sop'r'pi)=y^^]/2'^ 



(2) 



Of course, Equation [2] has no spatial reference, and therefore is not a complete representation of 
MFplp2; it gives only the frequencies of the richness at each level. The full MFplp2 is better 
thought of as a density where extremes of richness are rare, the richness distribution is unimodal 
at every scale, and near any cell with high richness there will be cells with low richness, but cells 
with high richness will tend to be near other cells with high richness. 



MFpl p2 is not very different from the 2-dimensional multifractal constructed by 



Stanley fc Meakin 



(|l988l ) : a much- cited example of the increasing use of multifractals in many fields including porous 



media 



Perrier et 



(|2006r ). At every second step in constructing MFplp2 one reaches a step in 
their construction. What is new here is firstly that pi and p2 have a natural interpretation in 
terms of species richness and secondly the explicit distribution of equation [2l 



A classical way to characterize a m ultifrac t al me asure is based on the calculation of the Renyi 
dimensions of a multifractal object. iHartd (j200l[ ) gives the definition as 



lim 



D{q) = { 



lim 



oq-l log2 r 



logs?- 



ifq^l 
if g = 1, 



(3) 



where r is a characteristic length, it takes iV identical boxes with length r to cover the object, and 
Pi is the proportion of the density in box i. 



Substituting the probabilities from Equation[2]into Equation[3]one obtains the following expression 
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Figure 3: Typical D{q) for v arious values of the variability scaling b. The flat line when 6 = 1 is 



the spatially invariant case of 



Harte et al 



(jl999l ) where pi — p2', SiS b grows the spatial variability 



becomes more pronounced and the curves depart more and more from a straight flat line. 



for the Renyi dimensions of MFplp2: 



2 , pI+pI 

l~q (pi +P2r 



if 9 7^1 



m = 



log2(Pi +P2) 



2 Pi logabi)^ -P2 log2(p2)^ 



(4) 



Pi + P2 



if g ^ 1. 

This in fact holds for finite fc, rather than merely in the limit of large k (see Appendix). 



Reparameterisation using a and h 



We now turn to an alternative and equivalent way to see MFplp2, via a — (pi + P2)/2 and 
b =pi/p2- 



Harte et al. 



( 19991 ) ■ as follows. 



Firstly, MFplp2 is a generalisation of the density constructed in 
When a = pi = P2, our construction reduces to theirs. Whenpi 7^ P2 > 0, we use the mean richness 
ratio a = {pi +p2)/2. Define mean richness Sk after k bisections as the sum of all the richnesses 
divided by 2'^. Clearly, the mean richness satisfies the recurrence Sk+i = {pi +p2)Sk/2 — aSk- 
This recurrence relation immediately implies the power l aw iS* = cA^, where z = logjl/a and 



Sq/Aq, exactly in the same way as in 



Harte et al. 



(jl999[ ). Hence the density of species 



richness they construct purely as a mean field model is exactly the trend of richness in MFplp2. 
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Secondly, the variability in S" at a given a depends entirely on b = pi/p2'. after k bisections, the 
ratio between the richest and poorest rectangles of size is b''. The variability increases by 

the factor b at every step, so b governs the scaling in variability. Interestingly, D{q) depends only 
on b: 



D{q) 



l0g2 (629 + 1 ) - log2 (6 +1)2? ^ ^ ^ 



1-q 

1 



In sum, a gives the scaling of the trend and b gives the scaling in the variability. Typical D{q) 
curves for MFplp2 are shown in Figure [3] 



Some properties of MFplp2 

1. Many statistics of MFplp2 have power laws. Above we derived S — cA'^, the usual interpre- 
tation of the species-area relationship. Similar expressions can be derived for the variance, 
minimum, mode and maximum of S{A). 

2. MFplp2 is exactly self-similar in the following sense. As k — > oo, every rectangle at level 

is a scaled copy of the initial rectangle, provided we consider richness relative to overall 
richness and ignore the placement of the two halves resulting from each bisection. 

3. MFplp2 has a characteristic triangular footprint on a log-log scatterplot, as in shown in 
Figure IH 

This illustrates graphically the way MFplp2 generalises the classical species-area power law: 
the triangular footprint collapses to the classic straight line when b — 1; the mean S of 
MFplp2 gives the classic straight line. 

In the example below, we show that it can match real data. 

4. MFplp2 predicts a positive correlation between endemism and richness. Let us suppose 
that Aq contains Eq endemics (species that do not occur outside Aq) and that endemics by 
themselves also divide according to pi and p2. Then the proportion of endemics that occur 
in the richer half with area Ao/2 is 1 —p2, because those species do not occur in the poorer 
half. That is, these species are endemics not only of the whole area Aq but also of the richer 
half. We define ei = 1 — p2 (respectively 62 = 1— pi) as the proportion of endemics that are 
endemic to the richer (respectively poorer) half of a a bisection, and denote by Pk{E ~ n) 
the probability that a rectangle with area Ao/2'' contains n endemic species. We then obtain 
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Log 2 (Area) 



Figure 4: The triangular footprint of MFplp2 on a log-log scatterplot. The non-vertical edges 
have slopes log2(pi) and log2(p2), the horizontal spacing is 1 {— log2(2)) and the vertical spacing 
is log2(6). The straight line is the classic Arrhcnius power law with a = {pi +p2)/2 and goes 
through the mean 5^ at every A. 



the endemics distribution Pk{E — E^el "'e^) — (^)/2'^. Because ei > 62 when pi > p2, the 
richest rectangles at a given k are precisely those with the largest number of endemic species. 

Just as for richness, statistics for cndemism have power laws. For example, for mean en- 
demism the exponent is ze — log^fl/g p), where ap, = 1 — a; here again we recover a result 



first given by by Harte and co-workers 



Harte fc Kinzij (|l997t ) 



Illustration: Proteaceae richness in the Cape Floristic Region 



The Cape Floristic Region (CFR) of the south-wes t ern co rner of Africa is the smallest of the 



planet's six floristic regions recognized in 



TakhtaianI (|l986l l. among which it has the highest en- 



demism. Proteaceae are the signature plants of fynbos, a vegetation type almost entirely confined 
to the CFR. The CFR is made up of several biomes, one of which is dominated by fynbos and hence 
is named the fynbos biome. Data o n species richn ess of Proteaceae was kindly provided by Tony 



Rebelo of the Protea Atlas Project 



Rebeld (jlQQlh . These data are derived from approximately 



60 000 site record sheets. These were filled in as follows. Firstly, the centre of the site was recorded 
using the GPS, then a roughly circular area with this centre and roughly 500 metres in diameter 
was designated, and finally all Proteaceae species found there were recorded on a standard form 
(as well as some other data irrelevant here; some pre-existing data were also used). 



In this paper, we consider a rasterization of the CFR by lines of longitude and latitude 1' apart. 
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Figure 5: Maps of richness of Proteaceae species in Cape fynbos. The first shows the finest 
available scale in the data: each cell is 1' x 1'. In the second, each cell is 7' x 7'. Cells are coloured 
according to number of species they contain, with the same colourmap from poorest (blue) to 
richest (red) in both maps; non-fynbos cells are blank. 

The resulting cells are very nearly rectangular, approximately 1 km by 1.5 km. The fynbos biome 
is here defined by protea presence: for the purposes of this illustration it consists of those cells 
that have at least one occurence of a Proteacaeae species recorded in the Protea Atlas. We call 
them fynbos cells, of which there are approximately 9 000 in the atlas. We wish to study the SAR 
of Proteaceae in fynbos. In other words, we need data giving area A and richness S samples from 
the fynbos biome. The samples were constructed as follows. Each fynbos cell was taken as the 
centre of a sequence of nested cells of size d minutes by d minutes. The smallest cell is of course 
1' X 1', the next size is 3' x 3' and so on. The full set of d are 1, 3, 7, 15, 31, 63, 127, 255, 511, 
1023. The fynbos area A inside a given cell is then simply the number of 1' x 1' fynbos cells it 
contains, and the count S is from the full list of species of these combined cells. Thus each 1' x 1' 
fynbos cell is at the centre of a nested sequence of fynbos areas, each with its own {A, S) values. 
At the largest scale, each cell contained all 9 000 fynbos cells. See Figure [5] for maps of richness 
using cells at two different scales. 

Figure El shows that the {A, S) data described above have a roughly triangular log-log footprint, 
and some values of pi,p2 that allow one to superimpose an MFplp2 footprint that fits it quite 
well. Of course, MFplp2 has only a small number of discrete values for (A, S) pairs, whereas 
the real data appear quite densely packed. Thus MFplp2 matches the pattern of the triangular 
footprint, but not all its details. It succeeds quite well in characterising the trends of both the 
mean richness and the variation in richness over many orders of magnitude of area. Using an 
algorithm for estimating Renyi dimensions from data (Perrier, E. and Laurie, H., unpublished 
work), the figure also shows that the data have a spectrum of Renyi dimensions very similar to 
that of an MFplp2 model. 



11 




Figure 6: Example of application of MFplp2: scaling of mean and variability of Proteaceae richness 
in Cape fynbos. (a) Triangular footprint of species richness, with an MFplp2 triangular footprint 
superimposed. Area is in number 1' x 1' cells; richness is number of species. Axes are scaled 
logarithmically. The blue dots are {A, S) pairs from the Protea Atlas, while the red dots are the 
triangular footprint of MFplp2 (with pi — 0.83, p2 — 0.61; equivalently a = 0.72, b — 1.36). (b) 
Spectra of Renyi dimensions: species richness of Cape Proteaceae (blue dots) and an MFplp2 
spectrum (red dots, pi — 0.83, p2 — 0.61; equivalently b = 1.36). 

DISCUSSION 



The variability of species richness at given area has been neglected in studies of SARs. We believe 
the reasons for this no longer hold. It is obviously very onerous to collect and process sufficient 
data to establish not only the trend but also the variability in S" at a given A. However, over 
time large data sets have been accumulated and with the developments in computing they can 
be handled. As t o theory, the construction of species density via an algorithm (pioneered by 



Harte et al 



19991) avoids two limitations of the derivation of SARs from SADs by probabilistic 
arguments. Firstly, no assumption is needed about spatial autocorrelation. As pointed out above, 
non-zero spatial aut ocorrelation emerges automa tically, and if need be additional assumptions are 



easily added (as in 



Martin fc GoldenfeL 



2 



20061 ). Secondly, algorithmic simulation allows direct 



modelling of processes at particular scales. By contrast, the allocation pro blem for the scale of 



19931 ): that is, from data 



spatial variation of smoothly varying densities is unsolvable (jCressie . 
alone one cannot identify the exact contribution of variation at each scale to the overall variation 
if the density varies smoothly. 



Let us note three recent papers that develop Hartean theory of spatial structure in SARs in 



ways quite different from MFplp2. All three embrace the SAD. 



Harte et al 



( 2005f ) introduces 



the concept of spatial abundance distribution, by which the overall SAD is distributed across 
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the cells created in a bisection cascade. Here, besides the parameter a of IHarte et ali (|1999l ). 
the abundances uq of each species are used. However, the traditional derivations do not apply, 
because of the spatial autocorrelation necessarily introduced by the bisection cascade and the 
parameter a. Instead, the analysis follows from a further assumption c alled HEAP, which we 

Martin fc Goldenfe J ( 2006 ): they also start 



do not treat here. Something similar is done by 
with the SAD as given and allocate the individua l s via a bisection cascade. However, instead of 
the parameter a and HEAP, iMartin fc Goldenfeldl (|2006r ) assume that individuals in each species 



are clustered. The generality of this assumption precludes analytical results, but their computer 
experiments are extensi ve and show c o nvinc ingly that realistic SARs do result from a variety of 
clustering assumptions. 



Conlisk et al. 



(|2007l ) is rather more preliminary. Here, the SAD is again 
the starting point for an allocation of species to a cascade of smaller and smaller regions. However, 
the geometry is generalised from bisection to n-section, whereby a given rectangle is divided into 
n >— 2 subregions at each step. The individuals in a parent region are distributed one by one to 
the subregions at each dissection step by a probabilistic rule. The authors show that this leads to 
coherent and reliable simulations, because for a given setting of the parameters, the final state of 
the system is independent of the exact sequence by which it was reached. They also study some 
statistical issues in parameter estimation with reference to three data sets. 

These papers differ from MFplp2 in three ways. Firstly, they treat spatial variation indirectly. 
Secondly, analysis of spatial variation is not pursued. Thirdly, they need at least as many pa- 
rameters as there are species. By contrast, spatial variation is directly modelled and quantified in 
MFplp2 using only two parameters. Thre usual trade-offs apply: loosely speaking, models with 
more parameters fit the data more closely, with less confidence in estimates of parameter values 
and more sensitive predictions. 

A particular strength of MFplp2 is that the parameter b = pl/p2 directly quantifies spatial 
variability across scales. To our knowledge, it is the first SAR model that quantifies variability 
independently of mean richness. We suggest that b may be a very useful summary statistic. 



The most obvious prediction of MFplp2, as clearly demonstrated by the log-log footprint, is that 
variability in S increases as A shrinks. This is in marked contrast with the one published resul t 



198l|) 



on spatial variability the we are aware of, namely the random placement model (jColeman 
where Var(5) goes to zero with A. The Proteaceae data shown above indicate that in at least 
some cases the prediction of MFplp2 is correct. 



13 



The ability directly to evaluate b as well as a from data sho uld allow comparis ons of variability in 



19951). It would be very 



S along the same lines previously used for the mean of S (jRosenzweig . 
interesting if situations are found where the values of a are similar while the values of b differ, or 
vice versa. It would be even more interesting if a and b were strongly correlated. 



We turn now to the objection (jMadduxl . 



20041 ) that bisection cascades are fatally impoverished: 
they furnish far too few values of A and far too few instances of domains with a given A. It is true 
that only k + 1 values for A occur from Aq to A^, and also true that only a few of the possible 
rectangles that can be formed from the recangles of size A^ that actually occur in the cascade. But 
the point is to describe the SAR, not to recreate it in detail by generating its value for all possible 
geometries at every A in a continuous range. That is, the aim is to summarise the behaviour of 
S versus A. In the case of Cape Proteaceae (see above), the log- log footprints of data and model 
quite clearly agree as to the trend of both the mean and of the variability. Here, MFplp2 appears 
to do a good descriptive job. The spectra of Renyi dimensions also match acceptably. 

The Cape Proteaceae example also shows that real data have quirks that a simple model could 
not possibly catch. Generalisations and extensions are necessary. We note that two possibilities 
remain quite simple. The first is that, instead of being constant, pi and/or p2 could be functions of 
A, thus introducing variation across scales. This might involve only a few extra parameters. The 
SAR would no longer be a power law, but the resulting density would still be multifractal (albeit 
not self-similar even in the weak sense used above) . The second is that the shape of Aq and indeed 
of all the Ai need not be rectangular: any shape, no matter how dissected or disconnected, can be 
bisected in many ways. Thus bisection cascades can be constructed for the typically complicated 
domains on which ecology is studied, and in many ways. 

Some immediate research tasks concerning MFplp2 are obvious. Application should be tested by 
estimating a and b for more sets of real data and studying the reliability of such estimates. We 
already have indications for several data sets that reliable estimates are possible (this is work in 
progress). The results of bisection cascades using various bisectors and rasterisations should be 
tested for consistency. There are also interesting open theoretical questions. One is the lack of 
a rigorous basis for inferences concerning SARs, because the correct probability distribution of 
domains with area A is unknown. The distribution cannot be based for all sets with area A, as 
the majority of such sets are ecologically meaningless. On the other hand, it cannot be based on 
just the convex sets, nor on the connected sets. Without a set of A to sample from, one cannot 
rigorously infer reliability, or indeed make any other inference. Note that this difficulty parallels 
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the problem of rigorous estima tion of the actual ( as opposed to observed) number of species, which 
is still the subject of research (jHaas et al\ . \200di) . 
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APPENDIX: Derivation of formula for D{q) 

First consider the case q ^ 1- By definition, D{q) — log^'/ ■ We need the pi, which are 

proportions of the total density. From equation [5] we know there are rectangles with richness 
Pi~''P2^o after k bisections. The sum of all richness is {pi +P2)'^So. In other words, 



— ^ — ( p\ ^p\ y _ (pI+pI)'' 



,_i ,_n \ (^'l +P2)'' / {Pl +P2Y'' 

I — 1 J — u \ / 



Each rectangle has an area proportional to 2 so one may use r = \Jl~^ as the characteristic 
length. It follows that /c — + cx3 as r — > 0. Hence we have 

k^ooq-1 log2 2-''-/2 fe^oog-1 -k/2 



We see that k cancels, so there is no need to take the limit for large k, and the first expression in 
equation [3 drops out. The second expression follows from the first by L'Hospital's rule. 



17 



